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ABSTRACT 

Observations of black hole binaries (BHBs) have established a rich phenomenology of X-ray 
states. The soft states range from the low variability, accretion disc dominated thermal (TD) 
state to the higher variability, non-thermal steep power law state (SPL). The disc component in 
all states is typically modeled with standard thin disc accretion theory. However, this theory is 
inconsistent with optical/UV spectral, variability, and gravitational microlensing observations 
of active galactic nuclei (AGNs), the supermassive analogs of BHBs. An inhomogeneous disc 
(ID) model with large (~ 0.4 dex) temperature fluctuations in each radial annulus can qual- 
itatively explain all of these AGN observations. The inhomogeneity may be a consequence 
of instabilities in radiation-dominated discs, and therefore may be present in BHBs as well. 
We show that ID models can explain many features of the TD and SPL states of BHBs. The 
observed relationships between spectral hardness, disc fraction, and rms variability amplitude 
in BHBs are reproduced with temperature fluctuations similar to those inferred in AGNs, sug- 
gesting a unified picture of luminous accretion discs across orders of magnitude in black hole 
mass. This picture can be tested with spectral fitting of ID models, X-ray polarization observa- 
tions, and MHD simulations of radiation-dominated discs. If BHB accretion discs are indeed 
inhomogeneous, only the most disc -dominated states (disc fraction > 0.95) can be used to 
robustly infer black hole spin using current continuum fitting methods. 
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1 INTRODUCTION 

Spectral and timing observations of black hole X-ray binaries 
(BHBs) in the energy range ~ 1 — 100 keV have established 
characteristic outburst "states" common to many sources (e.g., 
iRemillard & McClintockl l200d hereafter RM06). The "thermal" 
(TD) state has little variability, and a large fraction of the integrated 
flux is well described by a standard th i n black hole accretion dis c 
model (NT. Ishakura & Sunvaevlll973l ; iNovikov & Thornelll973l) . 
The "low/hard" state is characterized by a hard power law spectrum 
with a weak, soft disc component. It is often modeled as a thin disc 
truncated well outside of the black hole marginally stable orbit. In- 
side of this location, the accretion flo w is assumed to be geometri- 
cally thick and radiatively inefficient dEsin. McClintock & Naravan] 
fl99% . The "steep power law" state (SPL) is characterized by a soft 
(photo n index T > 2.4 ) power law extending unbroken to < 1 
MeV dGrove etaL 1998). Its physical origin remains uncertain. 

This common interpretation of the X-ray states relies on the 
validity of the NT model. However, this model is theoretically in- 
consistent in the TD and SPL states when the disc extends down to 
the marginally stable orbit of the black hole. In these luminous in- 
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ner regions, radiation pressure provides the ve rtical support against 
gravi ty, and the NT model is both thermally jshakura & Sunvaevl 
1 19761) and viscously dLightman & Eardlevll 19741) unstable. Global 
MHD simulations of thin discs dPenna et alll201(il ; ISorathia et alj 
I2012L and references therein), which include the physics of an- 
gula r momentum transport via the magnetorotational instabil- 
ity dBalbus & Hawley 119911 ). do not include radiation pressure 
and therefore cannot test the validity of the NT model. Local 
MHD simulations find tha t small patches of radiation-d ominated 
discs are thermally stable dHirose. Krolik & Blaesll2009l) . but sug- 
gest that the viscous (Lightm an-Eardley) instability may operate 
dHirose. Blaes & Krolikll2009h . This instability has also been pro- 
posed as a model for spectral state transitions in GRS 1915+105 
dBelloni et aHll997l ; lNeilsen. Remillard & Ledl2012l) . 

The NT model is disfavored by optical/UV observations of ac- 
tive galactic nuclei (AGNs) with similar L/L c< jd to the soft, lumi- 
nous states of BHBs. The m odel underproduces the observed UV 
emission dZheng et al.1 1997) and requires a relativistic mechanism 
to explain the simulta neous variability observed at well separated 
optical wavelengths dKrolik etalj|l99ll) . Recent microlensing ob- 
servations find that quasar accreti on discs are a factor of ~ 4 larger 
than predicted by t he NT model (jjimenez- Vicente et alj|2012l . and 
references therein). [Dexter & Agol ( 2011 , hereafter DAI 1) showed 
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that a disc with large, local temperature fluctuations (~ 0.4 dex) 
can explain all of these observations. The fluctuations could be 
driven by instabilities in radiation pressure dominated discs, in 
which case they would operate in BHBs as well. 

In this Letter, we propose a novel interpretation of the soft, 
luminous TD to SPL states of BHBs in terms of this inhomoge- 
neous accretion disc model (ID, summarized in Sj2j- We show how 
interpreting ID spectra with the standard NT model can explain the 
X-ray spectral and variability properties of soft BHBs ( |j3j> . In [j4] 
we assess the requirements for the ID model to fit observed BHB 
X-ray spectra. Inhomogeneity also may have important implica- 
tions (!|5} for the continuum fitting m ethod for measuring black 
hole spin (CF, lMcClintocketai]201ll. and references therein). We 
discuss the prospects for testing BHB accretion disc inhomogeneity 
observationally and theoretically in !j6] 



2 INHOMOGENEOUS ACCRETION DISCS 

DA 11 proposed a toy model of an inhomogeneous accretion disc. 
The disc is taken to be optically thick everywhere so that its flux 
is specified completely by the effective temperature, T(r,<f>,i). 
The radial temperature profile is set by assuming that the model 
reproduces the NT surface brightness when averaged over az- 
imuth and time. Retaining the NT model on average is moti- 
vated theoretically by conserva tion laws and numerical simulations 
I IHirose. Blaes & Krolikll2009f), and is consisten t with quasar mi- 
crolensing observations (Eigenb rod et al. 1 12008|) . DA11 broke the 
disc into independently fluctuating zones spaced evenly in log r and 
4>, with n total zones per octave in radius. Motivated by results from 
quasar variability (e .g. iKellv. Bechtold & Siem iginowska 1 20091 ; 
iMacLeodetail boiO). each zone performs a "damped random 
walk" in log T around the local thin disc value with amplitude <tt 
in dex and timescale tdrw Both of these quantities are assumed 
constant in radius, and the results are insensitive to tdrw- The re- 
sulting power spectra are broadly consistent with t hose observed in 
BHBs I Kell y. Sobolewska & Siemiginowskall201 ll) . but do not in- 
clude the high frequency quasi-periodic oscillations often seen in 
the SPL. 

For a fixed number of zones, the rms variability increases with 
increasing err. We fix n — 1000 (see § [3j>, in which case the rms is 
well fit by the empirical formula (cf. equation 1 of DAI 1), 



0.02 ^/exp (16a 



(1) 



In the limit n — > oo, this time-dependent model leads to a station- 
ary flux profile depending only on r and x = In (10) or: 
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dw 



-[lnw+x 2 ] 2 /x 2 



(2) 



where F v is the flux at frequency v and z = hv/kfT. The ef- 
fective temperature in each annulus is integrated over a log-normal 
distribution, whose width is In (10)ctt / 'V2. The local spectrum is 
assumed to be a colour-co rrected blackbody with / = 1.7 (e.g., 
IShimura & Takaharalll995l) . 

In the following sections, this stationary log-normal model is 
used for spectral calculations and the rms variability is estimated 
from equation ([]}. Using the log-normal model is valid for n > 30. 
Except in the brief discussion of X-ray polarization in fj(>] all rel- 
ativistic and viewing geometry effects on the observed spectra are 
ignored. Sample ID spectra are shown in the top panel of Figure[T] 
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Figure 1. Top: X-ray spectra from the inhomogeneous disc (ID) model with 
a = 0.5. The spectra become harder and broader with increasing or. Bot- 
tom: X-ray spectra from the NT (dashed) and NT(X>SIMPL (solid) models. 
A broader ID spectrum (dash-dotted line) with a power law tail (open dia- 
monds) can appear nearly identical to the NTXSIMPL spectrum. The power 
law components are shown as dotted lines. The black hole mass is IOMq, 
and the SIMPL parameters are f S c = 0.2, T = 2.7. 



3 THERMAL AND STEEP POWER LAW STATES 

RM06 define three states of BHBs in terms of their spectral and 
variability characteristics. The hardness ratio, the ratio of integrated 
flux from 8.6 - 18.0 keV to that from 5.0 - 8.6 keV, and the rms 
variability are model-independent. The disc fraction is the ratio of 
the integrated disc to total flux between 2 — 20 keV. The disc flux 
is found by fitting spectra with the NT model, as well as a power 
law tail component with a corresponding photon index T. 

Qualitatively, the ID models reproduce the observed trends in 
the TD to SPL states. For small or < 0.1, the spectra are nearly 
identical to NT spectra. They are soft, disc-dominated, and exhibit 
little variability (TD state). As <tt increases, the intrinsic disc spec- 
trum becomes broader and harder (Figure [I]. Fitting with a thin disc 
then requires an increasing degree of Comptonization, decreasing 
the disc fraction. The rms variability also increases with increas- 
ing <7t. For large or > 0.3, the ID spectrum extends to energies 
> 20 keV and bears little resemblance to an NT spectrum, so that 
the disc fraction is small. The rms variability has also increased 
significantly. This is similar to the observed SPL state. 

We demonstrate this quantitatively by treating the ID spectra 
as the data and least-squares fitting them with NT models along 
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Figure 2. Inferred spectral and variability properties of ID models with 
a = 0.5. Quantities are defined and plotted as in panels e), f), and g) of 
Figs. 4-9 of RM06. Each point is the result of fitting an ID spectrum with 
parameters in the range cry = 0.02 — 0.5 and rh = 0.8, 1.2 with the 
NTgDSIMPL model. The points are displayed according to their X-ray spec- 
tral state: thermal (red Xs), intermediate (open circles), and steep power law 
(green triangles). 
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Figure 3. Top: Black hole spin inferred from fitting a = ID spectra with 
the NT0SIMPL model as a function of disc fraction. If the ID model cor- 
rectly describes BHBs, the NT interpretation of the data can significantly 
overestimate the true spin for disc fractions < 0.95. The magnitude of the 
effect is sensitive to the properties of the high energy spectral tail. Bottom: 
<tt vs. disc fraction for the same fits. The disc fraction rapidly decreases 
for 0.15 < <Tt < 0.35, where the ID spectra become significantly broader 
than their NT counterparts. The ID spectra either have no high energy tail 
(blue filled circles), a power law matched on arbitrarily above 5 keV (green 
filled squares), or a power law tail from SIMPL (red open circles) with 
fsc = 0.03, r = 2.5. The error bars are the ranges inferred from models 
with varying rh at fixed <tt ■ 



with the SIMPL pr escription for Comptonization (NT®SIMPL, 
ISteiner et al.ll2009q) . With the simulated data and the parameters 
of the NT fits, we calculate the required RM06 parameters: hard- 
ness ratio and rms variability from the data (ID spectra), and the 
disc fraction and photon index V from the fits. The parameters of 
the ID (NT) spectra are a, rh, and <tt (a, rh, fsc, and T). For a 
given "source," the black hole spin is fixed, while rh and or vary. 
All of the NT parameters are allowed to vary. The normalized ac- 
cretion rate is rh = Mc 2 / L e dd- 

An example is shown in Figure [2] Each point represents a 
value of rh and <tt from an ID spectrum with a — 0, rh — 0.8, 
1.2, and a T = 0.02 - 0.5, in 0.02 intervals. The number of ID 
zones, n, is set by requiring that the rms variability not exceed 15% 
in the SPL state for a range of ID spin values. The inferred value 
n — 1000 is consistent with that required to explain optical/UV 
observations of AGNs (Figure 5 of DAI 1). The maximum value of 
or that leads to an SPL spectrum is ~ 0.5, identical to the max- 
imum allowed by AGN observations. Using a smaller value of n 
would restrict the range of allowed ot (rms < 15%), preventing 
the model from describing the lowest disc fraction SPL states. 

The results in Figure [2] are strikingly similar to the observed 
behavior of BHBs (e.g., compare with Figure 4 of RM06). The rela- 
tionship between disc fraction and hardness ratio is recovered with 
no free parameters. One point on the rms variability vs. hardness ra- 



tio relation is fixed by the choice of n, but the predicted trend is also 
in agreement with observations. Different choices for ID spin, or, 
and rh lead to variations of the same trends, and can explain much 
of the diversity observed in soft BHBs. For example, the presence 
of distinct tracks in plots of disc fraction vs. hardness ratio (e.g., 
Figure 7 of RM06) can be explained as variations in <tt at different 
values of rh. The inclusion of changes in n or possible correlations 
between n, gt, and rh could cause additional scatter in Figure [2] 



4 POWER LAW TAIL MODELS 

In the previous Section, the ID spectra were used as data with- 
out any power law spectral tail. This reduces the number of free 
parameters, but leads to steeper inferred photon indices than ob- 
serve d. The origin of this high energy tail emission is uncertain 
(e.g.. iDone. Gierlinski & Kubotd 120071 ). It can be modeled as a 
power law, often extending over all X-ray energies. However, at 
energies near the transition between disc and power law compo- 
nents, emission typically attributed to the power law could be due 
to a broader disc spectrum like that produced by the ID models. 
The high energy tail can also be modeled w ith physical 



Titarch uk 1994; Esin. McClintock & Naravanlll997h or empirical 



Steincr et al. 2009b) Comptonization models. In both cases, the 



disc spectrum (or its Wien tail) is used as the seed photon spectrum, 
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implicitly assuming that the Comptonizing corona extends uni- 
formly over the disc. This picture is strongly disfavored in AGNs by 
X-ray microlensing observations, where the c orona is found to be 
much more com pact than the optical/UV disc dChartas et a0l2009l ; 
iDai et alj20ich . Steep radial emissivity profiles found from model- 
ing X-ray reflection spectra in BHBs and AGNs (Miniutti & Fabian 
2004; ISchmoll et al.l i2009) also favor a compact corona. There is 
no theoretical reason to expect large differences in coronal struc- 
ture between AGNs and BHBs. All models for the high energy tail, 
then, introduce considerable uncertainty into BHB spectral fitting. 

This uncertainty may allow the ID model of the SPL state 
to fit the data as well as the significantly Comptonized NT one. 
An example is shown in the bottom panel of Figure [T] The same 
spectrum (solid line) is produced by a significantly Comptonized 
NT®SIMPL, and a strongly inhomogeneous disc supplemented by 
a power law with a low energy cutoff at 5 keV (diamonds). This 
power law could mimic Comptonization from a compact corona. 
Recent radiative transfer calculations from relativistic MHD sim- 
ulations have also found that emission from the plunging region, 
neglected in the NT model, can lead to a power law spectral tail 
dZhu et al.ll2012l) . Whatever the emission mechanism or geometry, 
the ID model for BHBs still requires that the normalization of the 
high energy tail increases with decreasing disc fraction (increas- 
ing <tt). ID spectra alone cannot explain the strong, non-thermal 
emission extending to 1 MeV in some SPL states. 



5 IMPLICATIONS FOR BLACK HOLE SPIN 

MEASUREMENTS VIA CONTINUUM FITTING 

One promising method for measuring black hole spin is to fit disc- 
dominated BHB spectra with the NT®SIMPL model. The hardness 
(or characteristic temperature) of the disc spectrum depends on its 
inner edge, which is assumed to coincide with the marginally sta- 
ble orbit of the black hole. This allows a measurement of the spin, 
provided that the mass, inclination, and distance are accurately de- 
termined. This tec hnique has been used to estimate the spins of 
several BHBs (e.g.. lMcClintock et alj20TTh . 

In the ID model, the characteristic temperature increases with 
both increasing spin and ctt. Therefore, if BHB accretion discs are 
inhomogeneous, the use of the NT model could introduce system- 
atic errors in CF spin measurements. From the fitting procedure 
used in Sj3] the magnitude of these errors can be estimated as a 
function of the inferred disc fraction. The results are shown in the 
top panel of Figure|3]for an ID model with a — 0. The range of or 
used is the same as in Figure [2] and is consistent with that favored 
by optical/UV AGN observations (0.35 < a T < 0.50, DA11). 

The inferred spin using the NT model is always larger than 
the "true" ID spin, since the spectral peak moves to higher fre- 
quencies with increasing or. Quantitatively, the spin differences 
depend on the high energy tail model used. This is because the spin 
and the SIMPL scattering efficiency fsc both harden the spectrum. 
Without any high energy tail (blue filled circles), fsc is small for 
<tt < 0.15 and instead the NT spin increases to fit the intrinsi- 
cally harder ID spectra. A power law matched on smoothly to the 
ID spectrum above 5 keV (green filled squares) leads to a higher 
power law normalization with increasing ctt. In that case, fsc must 
increase to fit the tail, and much smaller changes in spin can ex- 
plain the harder ID spectra. The disc fraction decreases sharply for 
°t > 0.15 as the ID spectra become broad, requiring significant 
Comptonization of the NT spectra in the fits (bottom panel of Fig- 
ure [3}. This leads to the sharp transition in inferred spin below disc 



fraction ~ 0.93 — 0.97, where the transition value depends on the 
high energy tail model. 



6 DISCUSSION 

X-ray spectral and timing observations of BHBs have established 
a set of states common to many objects. The physical origin of the 
SPL state and the cause of state transitions remain poorly under- 
stood. We have shown that the spectral and rms variability proper- 
ties of soft, luminous BHBs (TD to SPL states) can be explained by 
varying one parameter: the amplitude of local temperature fluctua- 
tions in an inhomogeneous accretion disc. Comparable temperature 
fluctuations in quasar accretion discs can resolve a number of im- 
portant observational puzzles in AGNs (DAI 1). The ID model then 
provides a unified picture of luminous black hole accretion discs 
across many orders of magnitude in black hole mass. 

The most likely physical processes for driving large temper- 
ature fluctuations are radiation pressure dominated disc instabili- 
ties. Local MHD simulations are consistent with the thin disc pre- 
diction for the relationship betwee n surface density and stress (o r 
temperature) when time-averaged dHirose. Blaes&Kroiikl [2009). 
The simplest interpretation is that the disc is viscously unstable. 
If so, it would likely produce temperature fluctuations qualita- 
tively similar to those in the ID model. This can be tested with 
radially-extended radiation MHD simulations. Alternatively, large 
fluctuations could also occ ur in magnetically-dominated discs (e.g., 
Begelman & Pringlel2007l) . where the disc radiation energy content 
is de-coupled from its stability. 

The ID model is insufficient to fit the X-ray spectra observed 
in both AGNs and BHBs. The broad ID spectrum produces some, 
but not all, of the emission typically attributed to the high energy 
tail, and some prescription for the remaining emission is still re- 
quired. This could be from Comptonization in a comp act corona 
or fro m hot gas inside the disc marginally stable orbit dZhu et alj 
l2012h . As with other models, the high energy tail normalization 
must increase with decreasing disc fraction in order to explain SPL 
spectra. Quantitative fitting of the ID spectral model to observa- 
tions will be carried out in future work, with various prescriptions 
for the high energy tail. Spectral models could use the colour- 
corrected blackbody used here, or could incorporate the results 
of sophisticated accretion disc atmospher e calculations (BHSPEC, 
iDavis et alJl2005l : lDavis & Hubemfcood) . 

If BHBs are indeed inhomogeneous, spin measurements from 
current CF methods using the NT model are potentially subject to 
large systematic uncertainties (spin differences Aa ~ 0.3 — 0.9 
for ID spin a = 0, Figure [3} . These can be much larger than cur- 
rent statistical uncertainties and systematic errors from ignoring 
emission from the plunging region (0.2 - 0.3. lKulkarnietalj201ll ; 
iNoble etalJkoilh . except in the most disc-dominated states (disc 
fraction > 0.95 or ot < 0.15). Therefore, if BHB discs are sig- 
nificantly inhomogeneous, spin measurements using CF methods 
should be restricted to the most disc-dominated states. 

The ID model also generically predicts that the inferred spin 
using NT m o dels should increas e with decreasing disc fraction. 
ISteiner et al.1 d2009al I2010L l201ll) have found that inferred spins 
in a few sources do not change significantly over a wide range of 
i/iedd or fsc- Since fsc is a proxy for disc fraction, this latter 
conclusion is in apparent confli ct with our Figure [3] (particularly 
Figure 1 of Stei ner et al. 2009a). However, the magnitude of the 
inferred spin differences when fitting NT models to ID spectra de- 
pends sensitively on the high energy tail model assumed. Further, 
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the inferred spin from fitting ID models is nearly constant outside of 
a critical range in disc fraction. Observations of IDs entirely above 
(below) this range in disc fraction could find a constant, correct (in- 
correct) spin value, although there is no evidence for sharp changes 
in spin from CF at disc fractions ~ 0.95. Fitting X-ray spectra with 
the ID model will allow estimates of a and ctt in a variety of states, 
and test the model. Requiring the spin to be independent of spectral 
state may constrain <tt and the degree of inhomogeneity in BHB 
accretion discs. 

X-ray polarization measurements may provide another test 
for inhomogeneous accretion discs in BH Bs. Assuming a semi - 
infinite scattering-dominated atmosphere (Chandrase kharl Il950h . 
we have calculated polarize d images and spectra via ray tracing 
using the cod es GEOKERR jDexter & Agol 2009]) and GRTRANS 
(Dexter 201 1]). We account for all r elativistic effects, including the 
rotation of the polarization vector dConnors. Stark & Piran|[l980l ; 
lAgolHl997h . With these assumptions, the time-averaged ID polar- 
ization is unchanged from the NT model: 1 — 5% peak polarization, 
increasing strongly from face-on t o edge-on viewing, and als o in- 
creasing at lower black hole spin (Schnittman & Krolik 2009). X- 
ray polarization may then provide an additional constraint on the 
BH spin, independent of the location of the spectral peak. In ad- 
dition, the polarization angle and degree vary by factors of a few 
on the fluctuation timescale. Both the degree of polarization and its 
time-variability increase with increasing or - These estimates ig- 
nore returning radiation, as well as coronal (|Schnittman & Krolikl 
and plunging region dAgol & Kro lik 2000) emission. These 
effects could cause additional differences between the polarization 
signatures of the ID and NT models. 

The ID model used here ass umes an optically thick disc. T he 
NT effective optical depth is (e.g., Ab ramowicz &Fragilell201ll) : 



TNT 



R 



/I 



93/32 



10M J 



(3) 



with G = c = 1 and ignoring the effects of general relativity and 
of the inner disc edge on the optical depth. If the disc is inhomo- 
geneous, portions will become optically thin. Assuming indepen- 
dent temperature and density fluctuations, we estimate the fraction 
of the disc with tid < 1 by integrating the log-normal distribu- 
tion, f{w), for tid(w) = w 2 - 3 tjs!t < 1. For a T < 0.1, > 99% 
of the disc is optically thick where the bulk of the luminosity is 
produced for L/L c( jd ^ 1- At larger <jt, the inner disc will be- 
come optically thin (> 50% by area at the peak emission radius 
for <tt > 0.35 and LjL B && = 1). If this em ission is similar to 
that from the plunging region dZhu et al .2012), it could lead to the 
observed increasing high energy tail emission with decreasing disc 
fraction. Comptonization from hot electrons in optically thin re- 
gions may als o effectively truncate th e disc spectrum as inferred in 
the SPL state foone & Kubot3l200d) . 



ACKNOWLEDGEMENTS 

JD thanks E. Agol, O. Blaes, J. McClintock, J. Steiner, and J. Tom- 
sick for useful discussions. EQ is supported in part by the David 
and Lucile Packard Foundation. 



REFERENCES 



Agol E., Krolik J. H., 2000, ApJ, 528, 161 
Balbus S. A., Hawley J. R, 1991, ApJ, 376, 214 
Begelman M. C, Pringle J. E., 2007, MNRAS, 375, 1070 
Belloni T., Mendez M., King A. R., van der Klis M., van Paradijs J., 1997, 
ApJ, 488, L109 

Chandrasekhar S., 1950, Radiative transfer. Oxford, Clarendon Press 
Chartas G., Kochanek C. S., Dai X., Poindexter S., Garmire G, 2009, ApJ, 
693, 174 

Connors P. A., Stark R. E, Piran T., 1980, ApJ, 235, 224 

Dai X., Kochanek C. S., Chartas G., Kozlowski S., Morgan C. W., Garmire 

G., Agol E., 2010, ApJ, 709, 278 
Davis S. W., Blaes O. M., Hubeny I., Turner N. J., 2005, ApJ, 621, 372 
Davis S. W., Hubeny I., 2006, ApJS, 164, 530 
Dexter J., 2011, PhD thesis, University of Washington 
Dexter J., Agol E., 2009, ApJ, 696, 1616 
— 2011, ApJ, 727, L24 

Done C, Gierlinski M., Kubota A., 2007, A&A Rev., 15, 1 
Done C, Kubota A., 2006, MNRAS, 371, 1216 

Eigenbrod A., Courbin E, Meylan G., Agol E., Anguita T., Schmidt R. W., 

Wambsganss J., 2008, A&A, 490, 933 
Esin A. A., McClintock J. E., Narayan R., 1997, ApJ, 489, 865 
Grove J. E., Johnson W. N., Kroeger R. A., McNaron-Brown K, Skibo 

J. G., Phlips B. E, 1998, ApJ, 500, 899 
Hirose S., Blaes O., Krolik J. H., 2009, ApJ, 704, 781 
Hirose S., Krolik J. H., Blaes O., 2009, ApJ, 691, 16 
Jimenez- Vicente J., Mediavilla E., Munoz J. A., Kochanek C. S., 2012, 

ApJ, 751, 106 

Kelly B. C, Bechtold J., Siemiginowska A., 2009, ApJ, 698, 895 
Kelly B. C, Sobolewska M., Siemiginowska A., 2011, ApJ, 730, 52 
Krolik J. H., Home K, Kallman T. R., Malkan M. A., Edelson R. A., Kriss 

G. A., 1991, ApJ, 371,541 
Kulkarni A. K. et al., 201 1, MNRAS, 620 
Lightman A. P., Eardley D. M., 1974, ApJ, 187, LI 
MacLeod C. L. et al., 2010, ApJ, 721, 1014 

McClintock J. E. et al., 201 1, Classical and Quantum Gravity, 28, 1 14009 
Miniutti G, Fabian A. C, 2004, MNRAS, 349, 1435 
Neilsen J., Remillard R. A., Lee J. C, 2012, ApJ, 750, 71 
Noble S. C, Krolik J. H., Schnittman J. D., Hawley J. E, 201 1, ApJ. 743. 
115 

Novikov I. D., Thome K. S., 1973, in Black Holes (Les Astres Occlus), 

New York: Gordon and Breach 
Penna R. E, McKinney J. C, Narayan R., Tchekhovskoy A., Shafee R., 

McClintock J. E., 2010, MNRAS, 408, 752 
Remillard R. A., McClintock J. E., 2006, ARA&A, 44, 49 
Schmoll S. et al., 2009, ApJ, 703, 2171 
Schnittman J. D., Krolik J. H., 2009, ApJ, 701, 1 175 
—,2010, ApJ, 712, 908 

Shakura N. I., Sunyaev R. A., 1973, A&A, 24, 337 

— , 1976, MNRAS, 175, 613 

Shimura T., Takahara E, 1995, ApJ, 445, 780 

Sorathia K. A., Reynolds C. S., Stone J. M., Beckwith K, 2012, ApJ, 749, 
189 

Steiner J. E, McClintock J. E., Remillard R. A., Gou L., Yamada S., 

Narayan R., 2010, ApJ, 718, LI 17 
Steiner J. E, McClintock J. E., Remillard R. A., Narayan R., Gou L., 

2009a, ApJ, 701, L83 
Steiner J. E, Narayan R., McClintock J. E., Ebisawa K, 2009b, PASP, 

121, 1279 

Steiner J. F. et al., 201 1, MNRAS, 416, 941 
Titarchuk L., 1994, ApJ, 434, 570 

Zheng W., Kriss G. A., Telfer R. C, Grimes J. P., Davidsen A. F, 1997, 
ApJ, 475, 469 

Zhu Y, Davis S. W., Narayan R., Kulkarni A. K., Penna R. E, McClintock 
J. E., 2012, |arXiv:1202.1530l 



Abramowicz M. A., Fragile P. C, 2011, arXiv: 1104.5499 

Agol E., 1997, PhD thesis, University of California, Santa Barbara 



© 2012 RAS, MNRAS 000,[TH5] 



